UNIVERSITY OF LONDON 
INSTITUTE OF ARCHAEOLOGY 


DENUDATION CHRONOLOGY 
_THE DYNAMICS OF RIVER ACTION 


ony: | 
DAY KIMBALL 
_ With an Introduction 
by : 
F. E. ZEUNER 


; 5 “PUBLISHED AS OCCASIONAL PAPER No. 8 
4 ar THE INSTITUTE, INNER CIRCLE, REGENT’S PARK, LONDON, N.W.1 


PRICE TWO SHILLINGS AND SIXPENCE 


1948 


UNIVERSITY OF LONDON 
INSTITUTE OF ARCHAEOLOGY 


DENUDATION CHRONOLOGY 
THE DYNAMICS OF RIVER ACTION 


BY 
DAY KIMBALL 


With an Introduction 
by 
F. E. ZEUNER 


PUBLISHED AS OCCASIONAL PAPER No. 8 
AT THE INSTITUTE, INNER CIRCLE, REGENT’S PARK, LONDON, N.W.1 
1948 


~ CONTENTS 


INTRODUCTION, BY F. E. ZEUNER 
Part I. UNGRADED STREAMS 


Part II. GRADED STREAMS 


PAGE 


DENUDATION CHRONOLOGY 
THE DYNAMICS OF RIVER ACTION 
By Day KIMBALL 


INTRODUCTION 


PrEuisToric Chronology is closely linked with the chronology of rivers. The 
Thames and the Somme are outstanding examples of rivers whose “ terraces ”’ 
are used in assigning relative dates to palaeolithic industries, and the Ilm-Saale 
system in central Germany is an example of a river system which has been 
used in the reconstruction of a detailed climatic chronology for the Pleistocene. 
It is curious that, although the importance of the activities of rivers in 
Pleistocene chronology (and in its application to prehistory) is readily admitted, 
very few authors have made any effort to discover what rivers do when certain 
environmental changes occur. Unfounded assumptions are often met with 
(such as the conception of “ overloaded ” rivers which Mr. Kimball explodes 
in the present paper), and others with only slight foundations in fact (such as 
the conception of aggradation in arid climates), which have not been studied 
in detail. : 

Clearly, if we want to use for Pleistocene and Palaeolithic chronology the 
periods indicated by the alternation of erosion and aggradation, or of intense 
erosion and erosion reduced almost to zero, as they are exhibited by many 
rivers, it is imperative to analyse the behaviour of flowing water. This is a 
subject which has for a long time occupied the water-engineers, but their 
interest has been not so much in natural rivers as in rivers that have been 
interfered with by man, and their main question, how one must interfere still 
further in order to make the river comply with man’s requirements. 
Consequently, though such investigations have produced many valuable 
results, they do not answer the questions which we, as geologists or pre- 
historians, are asking. ‘These questions have not been tackled in more than 
a passing manner since the days of the great American geologist, G. K. Gilbert, 
who published his Transportation of Debris by Running Water in 1914. His 
results are still regarded by many as the final word in these matters. Now, 
we owe a debt of gratitude to Mr. Day Kimball who has for some years studied 
with great care the relevant papers of earlier workers. In particular he has 
undertaken to interpret the results of water-engineers in terms of physical 
geology, by no means an easy task as the reader will notice when he continues 
to peruse this paper. We are indeed fortunate in having in Mr. Kimball, 
a member of the Institute of Archaeology.who is capable of coping with the 
exceedingly complex problems involved, and one who has spared himself no 
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efforts. to understand the mathematical methods of representation used by 
the various authors. 

The nature of the material under discussion makes the use of mathematical 
symbols inevitable. Mr. Kimball has developed ingenious methods of 
obtaining simple formulae for the relations of factors contributing to river 
action which have escaped earlier authors. His approach is new in several 
respects and the results, from the geological standpoint, gratifyingly simple. 
But as I fear that the majority of readers are liable to stop reading when 
encountering the first mathematical symbol, either from respect or disgust, 
it is perhaps useful at this point to say in a few words what the paper is about. 

The chief characters of a river are (1) its discharge, i.e. the quantity of 
water passing the point of observation, measured in cubic feet per second ; 
(2) the slope of the water-surface, usually referred to simply as slope or 
gradient of the river ; (3) the mean velocity of the flow ; (4) the shape, 1.e. the 
cross-section, of the bed ; and (5) the load carried by the water. Of these, 
(5), the load, is not treated in the present paper, for reasons given on p. 13. 
The others are defined on pp. 4-5, in particular the elements composing the 
cross-section. 

In Part I of this paper, Mr. Kimball develops formulae which enable one 
to calculate, very approximately, what happens to a river when the discharge 
is changed. In other words, if the climnate of an area changes, and the rainfall 
increases or decreases, what effect would this change have on the velocity 
and the shape of the bed? As a concrete example the Thames is taken and 
it is found that, if Atlantic times were to return to England, the immediate 
effect on the Thames would be to quadruple its discharge, to increase the 
velocity by more than one-half, to double the depth of its channel and to 
increase its width by a quarter. 

But the formulae developed in Part I only give us the zmmediate changes 
that would occur. The slope of the river, for instance, would not be affected 
straight away, and the river would be in an unbalanced condition in which it 
would try to adjust its bed and its longitudinal profile to the new conditions, 
by down-cutting and/or deposition of gravel. Such rivers are called ungraded. 

In Part II, therefore, the question is raised—and answered by approxi- 
mative formulae, what the ultimate result of such a change in annual 
precipitation would be, in other words what the slope and the shape of the 
channel, and the velocity will be when the river has remodelled its bed and 
thus has become graded. 

Archaeologists who are familiar with the numerous Palaeolithic sites 
found on the banks of ancient rivers and which are preserved in “ river 
terraces’ will appreciate the significance of the formulte developed by 
Mr. Kimball. They will help both in making more definite statements about 
the environment in which early man was living and also in placing terrace 
chronology on a firmer basis. 
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PARTS] 


UNGRADED STREAMS | 


Historica, Geology may be divided into two parts : Stratigraphy and its 
baby sister Denudation Chronology. The data ‘of the two differ somewhat. 
The first deals with what is there : the second with what isn’t. 

Not unnaturally, therefore, the curse of the business is a “ subjectivity ”’ 
In reconstruction far in excess of that which afflicts the older Stratigraphy. 
To this I am confident anyone will agree who has ever struggled to link up 
the denuded and dissected remnants of a flight of river terraces, and has 
struggled at the same time to retain some few shreds of “ objectivity ” in his 
reconstruction. 

Now Denudation Chronology is largely concerned with the. vestiges of 
past fluviatile processes. Here, surely, is a hint of one way to achieve greater 
objectivity ; which is to learn more than we geologists now know (and that 
should not be difficult !) about the mechanics of river action, about what an 
engineer calls the hydraulics of flow in an open channel. 

Clearly what we need is the answer to such questions as: Which are the 
fundamental variables ? Which of these are dependent and which indepen- 
dent? What factors, exactly, determine each of the variables, and in which 
sense, directly or inversely ? What is the order of magnitude, at least, of the 
effect on each variable of a given change in one or more of the factors that 
determine it? For guidance in such matters we must turn to the Engineers 
and Physicists. 

The purpose of the present paper is to attempt to answer these questions, 
by recording, in a form acceptable I hope to geologists, some few nuggets of 
truth that I acquired while trespassing shamelessly in the field of a sister 
science. ‘Thanks to a congenital inaptitude for mathematics, it cost me much 
painful slogging to disengage them from a ground-mass of less relevant truths. 
For this, however, I shall feel rewarded if they can be deemed a first step 
towards greater objectivity. Yet even if they can, that step is—like the 
housemaid’s baby—a very little one. 

First, then, as to those dependent or “internal” variables that are 
completely within the control of the river. They will reflect immediately any 
change in one of the “‘ external” variables" : 


W = the width of the water surface, in feet. 
D = the maximum depth of water, in feet. 
A = the cross-sectional area of the stream, in square feet. 


(1) Such, for example, as a change in the regional rainfall, and hence in the size of the river. 


UNGRADED STREAMS 


A 
== the mean depth = W 


D 

V = the mean velocity of flow, in feet per second. 
P= the “wetted? permmeter. = injfect. sioce.Mig.ar.) 
R 


= the “ hydraulic mean depth’, in feet = . 


c ol e 
Fic. 1. 
Diagrammatic cross-section of a river, to illustrate some external variables. 
D=h-d. ; 
® W=a-g. 
P=a-b-c-e-f-g. 


As drawn W=5D. 


Of these seven, the “‘ fundamental variables”? are V, R and P.* But 
fortunately, as geologists, we need hardly ever distinguish between P and W ; 
or D and R. Provided W exceeds 10 D®, P will exceed W by less than 5 per 
cent.) ; and since A/P = R, and A/W = D, D will exceed R by the same 
amount. ‘The present paper deals solely with streams in which W exceeds 
10 D. 

We may, therefore, in our quest for orders of magnitude, use the observable 
W in place of P. The relation of D to the observable D, is purely geometrical. 
If the channel is rectangular D = D; if triangular D = 3D. Im all other 
forms of channel the value of D lies between these two limiting values. Unless 
we know, or have assumed™, that our channel has one of the limiting forms, 


we shall unblushingly write the approximation: R = D = 2 D. 


The independent variables are : 


Q = the discharge, in cubic feet per sec. 
S = the slope of the water-surface, in feet per foot". 


(t) Lacey, 1935, p. 536. 


(2) In the Thames above Teddington, at bank-full stage, W=20 D, and P only exceeds W by 
1-6 per cent. ; 


(3) Lacey, 1930, p. 290. 
(4) cf. infra, pp. 8 and 9g. 


(5) We are primarily concerned with the river in its highwater stages. Hence S is, for us, the 


slope of the surface of the floodplain. It is assumed throughout that this surface and the mean 
bed-slope are sensibly parallel. 
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Here a distinction must be made. Our problem requires the proverbial 
two bites. The first question is: What will be the immediate consequences 
of a change in one of the independent variables? If we designate by Q.,, 
V,.. etc., the original values of the variables; and by QO); V,.. . ete., 
their values immediately after the change; then a change from Q, to Q, 
will cause no immediate change in S. S, will remain equal to S,. For a river 
cannot change its bed-gradient appreciably, in an hour, or even a day. 
S, therefore, is an independent variable, like Q ; whereas all seven variables 
in the first list are immediately affected by a change in S or Q," 

On the other hand, ultimately, a river “ grades ” its own slope. We may 
designate by S, the ‘“‘ graded” slope towards which a river is constantly 
tending and which (in theory at least) it will ultimately attain, if given 
sufficient time. Similarly Vg, We . . . etc., designate the ultimate ‘‘ graded ” 
values of V, W, etc. The point to observe is that, in calculating Vo, W. .. . 
_etc., S is no longer an independent variable but has become a dependent 
variable like V or W. 

A third factor, also, behaves in a similar fashion. The relations between 
W, D, D and A are purely geometric. They depend upon the proportions and 
shape of the channel. This (like S) is an independent variable as regards the 
immediate effect of a change. It takes time for a river to alter its channel. 
But the ultimate, “graded” shape is fashioned by the river itself, and, 
therefore, becomes a dependent variable. Here the point is that the initial 
geometric relations between W,, D,, D,, and A, may be regarded as constant in 
calculating V, or W, ; but may not be so regarded if we are seeking Vor We. 

The simplest, and therefore, for our purpose, the best of the various 
“ Flow Formulae” is Manning’s (Gibson, 1930, p. 288). It is usually 
written : 


Te40G Rests fr st fe ee: (1a) 


(t) Personally I find a further distinction helpful. Unlike a mathematician, perhaps simply 
because I am such a bad one, I refuse to regard even the most seductive differential equation as an 
ultimate : for, as Lord Rutherford used to Say, 1t ignores ‘‘ the real cause of the thing’ (cf. Eve’s 
Rutherford, pp. 208 and 221). To “understand” a natural process means (for me) to be able to 
construct a mental “ model ”’ of it and to state in words how it works. (I refer to observable, 
macroscopic processes.) Now the snag in picturing hydraulic processes is the deplorable way in 
which the dependent variables mutually determine one another. To insist, as I have done, that, 
physically, Q and S are arguments not functions, helps, but does not end one’s difficulties. V 
unquestionably (in part) determines D. But it is also, and equally, true that D (in part) determines 
V! The solution (for me) is to postulate, quite arbitrarily, that all geometrically determined effects 
are “instantaneous ’’, whereas all physically determined effects entail a slight time lag. By 
substituting two simple but successive processes for one complex process, we falsify nothing : for we 
come out, after two steps, at the same point that nature reached in one. The procedure is strictly 
analogous to that invaluable device of Physicists, the Resolution of Forces ; except that my 
“resolution ”? is temporal instead of spatial. It permits one to say (or think) : “To double Q will 
instantly double A and, assuming a rectangular channel, D and D also. Hence D, will=2 D,. Then 
the consequent increase in V will progressively decrease this (hypothetical) D,, until equilibrium is 
attained at Dz, the point where... . etc., etc.” My second distinction, therefore, is between 
“instantaneous ” geometrical and ““immediate ” but “ non-instantaneous ” physical, effects. 
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For normal rivers, like the Thames for instance, N has a value of about 
0°0225 3. Siving = 


V = 65 R' S? (approx.) .. ie re ut (tb) 
But for us the most useful form is : 
Wiroc Ry Stees 2 a c; a A: (ac) 


The only other equation we need") is that which expresses the fact that 
we are concerned with “ steady ’’, not “* fluctuating ”’ flow : 


OF VAT VERS wi D> a ie. af: (2) 

The first point to notice is that equation (1c) alone does not permit us to 
determine the amount by which V will be increased, if, for example, S be 
doubled. This must be stressed because the contrary is, apparently, often 
assumed, 

The physical significance of the formula would be better expressed by : 

e 3 
This emphasizes the fact that how much of the increase in S will be 


reflected by an increase in V, and how much by a decrease in R, is no concern 
of the Flow Formula. 


* The “allocation”? between V and R of the immediate consequences of 
a given change in 8 (or Q) is determined by the shape of the channel. 
Consider first a rectangular channel : 


oc S? 


Here W will remain constant. Hence, from (2), since P = W: 
OVE EB RecoaV hinace by oc 


Substituting in equation (Ic) : 
g 
Wee ROSE ee = Ss? 
“Vo Q'S 
Hence, if we designate by QO* and S* the factors by which we know (or 


assume) that Q, and S, were multiplied (to give Q, and S,) and by V* the 
factor by which V, will be multiplied (to give V,), we may write : 


V¥ = QO*SRR . ae is (3) 


Also, since R oc Vv: 
(1) Kimball, 1939, p. 144. 


(2) Thus we find in even so admirable a text book as that of Longwell, Knopf and Flint, end ed., 
1939, p. 61, the statement : “ If the gradient be increased four times, velocity is about doubled.” 
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UNGRADED STREAMS 


Se ones 
Q: Sto Sito 


D* =e ae ie a =f i a (4) 


; and we may write : 


eo 
Ow 


Finally, since W remains constant : 


Consider now a triangular channel : 
Here W will not remain constant but will vary directly as D. And 
D =D. Hence: D also varies directly as D. Hence: 
JP c= MV we ID ee ID) alte. 
Hence, from (2) : 
OF PReoe VV RRe oe V Re 
Q:? 
Substituting, as before, in (Ic) : 
V cc Sh Be NCE ORS BE NY we (OES Bate 
Vices Ota S ke a - = ig me (6) 


Q} a 
Also, since R oc, and D= 2D 


Vv 
ED Ne aT ee Ee 


Tees ae! a Ae aE is Meaiceos (97) 


OL: 


«Ors 
Also, since A o OF 
V 
Q Oe. 
Ao Ors: nC or : 
Ox 
Finally, since Po R: 
W=Pa es: aye 
O*s : 
NE rrr Fc or a bes me oe: (9) 


Now the two forms of channel just considered are the two extreme, or 
limiting, cases, All natural channels lie between the two. Here at last is 
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something to get our teeth into. For we can now say of any natural channel 
whether we know its shape or not: the immediate effect of multiplying 
Q, by Q* (or S, by S*) will be to multiply V, by a V*, the value of which will 
lie between the values given by equations (3) and (6). Further, since it is 
order of magnitude we are after, we may take the mean of the two values as 
the value of V*—approximately. Similarly the value of D* will be, 


approximately, the mean of the values given by equations (4) and (7) ; 
A* 7 
and W* will have the approximate value =;,, where A* is the mean of the 


De 
values given by equations (5) and (8). 
But we need not stop there. Our aim being what it is, we may go one step 
further. Instead of averaging the values given by our paired equations, 
we may “average ’’, as it were, the equations themselves. We then have, 


at last, what we really need: the four very approximate, but very simple 
equations : 


Vie Orr ot PRSOH Cr ius em heey i gen Ona) 
gy eee | | 
ears ms ous he ors iA oes (11) 
At oe 7 a « af a ee (12) 
Q* 
Ww* Cie °° es “i ie ie es (13) 


I shudder to think how a mathematician would characterize such 
concoctions. Yet I believe that, because of his very special needs and purposes, 
they give a Geologist something of considerable value”. Let me illustrate 
by a concrete example. 

The land snails of our Holocene deposits speak unambiguously. During 
Atlantic times the precipitation in the Thames basin must have been far 
greater than it is to-day. A conservative estimate would put the precipitation 
at Kew, then, as equal to that of Devonshire or Cornwall to-day. This 
means that it was about double what it is to-day. If we were to double 
the precipitation in the Thames watershed to-day, we should also 
approximately double the percentage of run off’). We should therefore 


(t) The cause of the mathematician’s sturdy refusal to supply us with such blunt but essential 
tools is presumably that he does not grasp the special needs of the Geologist and can not conceive 
that such shoddy and discreditable implements could be of use to anyone. It was this refusal, by 
Mathematician, Physicist and Engineer alike, to give us what we need, that drove me to the foolhardy 
attempt to steal the necessary material and fabricate the things, as best I could, myself. Why has 
mathematics no branch called ‘‘ Approximations ” lying somewhere between “‘ Applied Mathematics ”’ 
and ‘‘ Probability ” ? 


(2) The Weekly Weather Report, vol. LV (N.S.) (1939), pp. 46, 59, 60. 
(3) Thames Conservancy Statistics, London, 1935, pp. 16 and 17. 
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multiply Q, by about 4. .. Q* = 4. Hence, from equation (10), (11) 
and (13): 


V¥ = W4 = 1-59. 
IDE ces Gy) tee ty 
yA AT 20, 

We may, therefore, assert with some confidence that, if Atlantic times 
were to return suddenly to England, the immediate effect would be to quadruple 
the mean discharge of the Thames ; to increase its present mean velocity by © 
over 50 per cent. ; and, where the banks were high enough, to increasé its 
average mid-channel depth by 100 per cent., and its average width by 
25 per cent.” 2 

No one is more conscious than I of the many shortcomings and demerits 
of these results. One outstanding defect is obvious. All we have got, as yet, 
is the zmmediate effect of the change ; what would help us far more would be 
to know its ultimate effects ; for these would determine what the Thames was 
like (or, more accurately, what its graded reaches were like) towards the end 
of Atlantic times. But we have taken a step. That is something. And 
though still far short of our objective, what we have got has, in my eyes, one 
feature of rather exceptional merit. If the initial assumption was justified, 
then the results are as objective as those of the Physicist. True he desires and 
achieves extreme accuracy ; while our results give us only orders of magnitude. 
But then that form of knowledge is what, and all that, really concerns us as - 
Geologists. , 

I shall not attempt to carry the matter further in the present paper. 
I am confident that it can be done. Whether I can do it is another matter. 
But, at some future date, I propose to try my hand at it. Meanwhile it is 
no bad thing that the fundamental difference between the two problems 
should be emphasized by a break ; by a minor unconformity as it were. 

That there may be no possible ambiguity as to the limits of what I think 
I have accomplished in the present paper, I shall, in conclusion, emphasize 
three points. 

(1) Equations (3) through (13) do not purport to yield absolute values 
of so many feet or feet per second. They are designed to give the approximate 
factor by which the original value of the dependent variable will be multiplied 
as the result of a given change. Take for example equation (10). 


(t) Had we, instead of using equations (10), (11) and (13), used our paired equations, we should 
have got : 


From (3) and (6) : ven 57, 
; : -68 
From (4) and (7): D*¥=7 or —1-99. 
: 8 A® 9-56 
From (5) and (8) : ie eo Se oan giving for We Ds ; aga 28. 
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V*¥ = Q*! $**, What this asserts is that V,: V,77O7: QO.) Henceat OO, 
is a given number of times greater than Q, and we call this number Q*, 
then O,=Q* x Q,; and QO = O*F x OQ}. Substituting in the original 
: Ver OR OF rare — (j*t 

proportion : Wie le ie FOr sO) tec eee lirica hem actor 
by which V, will be multiplied is Q**. And this we have expressed by writing 
V* = Q*? §**, In the. case supposed 5S, remained constant., 5, = 5,. 
Hence, S* = 1. IfS had varied the same reasoning would apply to S*. 

(2) The equations deal solely with the immediate effects of a change in 
OQ or S or both. They do not purport to tell us anything about conditions after 
the river has begun to adjust itself to the initial change. If the reach in 
question was graded before the change, the process of adjustment will tend to 
restore it to a graded condition. If there is time enough, and, if no new 
factors intervene”, the reach will once more become graded. When, and if, 
this happens, a new set of equations will apply. They will yield values 
fundamentally different from the equations of the present paper. ‘The new 
equations are more difficult to derive, because 5 and “shape ”’ must now 
both be treated as dependent variables. But provided that Approximations 
are sufficient for our purpose, I am certain that such equations can be 
puCONnCOCLEC. 

(3) The equations contain no term that purports to reflect the amount 
of load carried. ‘This does not mean that the equations are restricted to an 
unloaded stream. Apart from exceptional cases”, all natural rivers carry a 
load ; and most carry the maximum load which they can transport. 

Shortly before the war the Thames Conservancy completed a project 
that practically doubled the cross-sectional area of the Thames for 11 or 
12 miles, from Shepperton to Teddington. It was necessary to calculate in 
advance the effect of thus doubling A upon the other variables—especially S. 
The formula used was Bazain’s. ‘This differs from, but is analagous to, the 
formula (equation 1a) from which I have derived all my approximate 
equations. 

The value selected for the constant in the Bazain formula would have 
been represented in Manning’s formula by: N = 0-0225. ‘The water level 
at Shepperton was, in fact, lowered several feet. ‘The actual amount differed 
from the calculated amount by from 2 to 3 inches. 

Someone had blundered (in calculating S,), . . . by almost 0-0000035 
units !'9) 

When the Thames reaches Teddington it is carrying a very substantial 


(t) Such, for example, as the retrogressive erosion of a knickpoint from some lower reach into the 
reach in question. 


(2) The Rhine from Lake Constance to the falls at Schaffhausen carries practically no load. 
(3) Naturally one must remember that, with a slope of one foot per mile, So: 0002 feet per foot 
approx. But putting it as I have above makes a better story of it. 
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load”). Equations (3) through (13), therefore, like their parent formula (1a), 
are most emphatically applicable to a loaded stream. What they do not do, 
is to enable one to treat “‘ Load”, also, as an independent variable and to 
calculate the approximate effect on the dependent variables of a given change 
init. In other words, it may be, though this is by no means certain, that our 
equations rest throughout on a “‘ concealed premise ” that : (Qathe. echarcecs 
(in Ibs. per cubic foot of water) and (b) the mean “ grade’ of the material 
transported, shall remain substantially constant. 

How the load of a stream can be expressed mathematically ; and how a 
change in load will effect the dependent variables is, at the moment, a highly 
controversial topic. Whether some day we can add to our equations a term 
representing load is not yet clear. If we can, the new term will have to les 
it would seem, a composite one®’, the two components reflecting in some form 
the two complementary aspects of ‘‘ load”, as a Geologist understands that 
term. For what does seem factually certain is that the mass transportable 
varies inversely as) the grade of the detritus. 


(*) See Thames Conservancy Statistics, p. 6. 
(2) i.e. of some such form as: (4a xb)x. 
(3) Though not, necessarily, proportionately to. 


8) 
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GRADED STREAMS 


WE now turn to a new and wholly different problem : If a stream, or a reach, 
be graded, what will be the wultzmate') effect upon the various dependent 
hydraulic variables”) of some such external change as a permanent increase 
in the mean annual rainfall of the catchment area ? 

If the new conditions persist substantially unchanged for a sufficient 
length of time, each of these variables will assume and maintain’) a new value. 
The river, or reach, will then once more be graded. What we are after is a 
petro paaors that will tell us (approximately) how these new and 
permanent" values are related to the values which the variables had before 
the change occurred. Will the slope, for instance, be twice what it was before, 
or half what it was before, and so on. : 

The equations will be similar to those derived in Part I, but the physical 
reality for which they stand is profoundly different in the two cases. Each 
of the equations of Part I yields the approximate value of an actual, immediate 
change that could be verified by appropriate measurements. Each of the 
new equations will yield an approximate value that is hypothetical and in 
a sense, ideal. It is the value which that particular variable will tend to 
approach more and more closely, decade by vanes and which (in theory at 
least) it will ultimately attain and maintain"), provided no further external 
change occurs in the meantime. Yet teguiles. their hypothetical and ideal 
character these ultimate “‘ graded” values are concepts which possess, I 
believe, greater utility for the geologist than do their more ‘“‘ real ’’ but more 
ephemeral counterparts. 

‘In the new equations the symbols will be similar to those employed in 
Part I, but with the subscript G added. Thus : 


Vee wineall evelocity 

V, = the initial mean velocity before the change. 

V, = the mean velocity immediately after the change. 
Vo = “the ultimate, eraded sincant velocity. 


V* = the factor by which V, was (or must be) multiplied to yield V,.. 
V*,, = the factor by which V, was (or must be) multiplied to yield V<. 


(1) See ‘p. 7. 
imocer pans: 
(3) Or perhaps more accurately : assume and thereafter oscillate about . . . 


(4) Permanent, that is, so long as the external factors remain substantially constant ; and ignoring 


the almost imperceptible, secular diminution in slopes that will accompany the pees ‘peneplanation, 
of the region as a whole. 


(5) See notes (3) and (4) supra. 
(6) See pp. 5 and 8. 


we 
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GRADED STREAMS 


My original intention had been to attempt, in Part II, to derive my own 


approximate equations for S*., V*o, D¥, . . . etc., by the same sort of method 
as that employed in Part I. 


It now seems to me preferable to adopt a different course ; I shall, 
instead, adapt to the purposes of Geology and translate into the notation 
already adopted in Part I, a few of the numerous equations derived by Colonel 
Gerald Lacey in two notable papers contributed by him to the Proceedings 
of the Institution of Civil Engineers. 

Despite their arresting character and the obvious value for Geology of 
much that they contain, it at first seemed to me, impossible to use these two 
papers directly’). One thing only stood in the way. ‘That “ thing”? was an 
objectionable small entity called by Lacey “f”. It appears as a variable 
in nearly all his primary equations. Now Lacey’s equations are all, 
ultimately, based upon empirical data. “‘f’’, therefore, has no a priori rational 
significance ‘“‘ by definition’). In some sense or other it represents the 
stream’s “‘ load ”*S), | , 

But precisely what it represents, physically, is by no means clear. Nor 
is there any agreement on this point among Lacey’s critics. In fact the only 
point upon which they appear to agree is that Lacey is wrong in insisting—as 
he does—that “‘f”’ is determined by, and solely by, the mean grade of the 
bed-load). 


It was this unfortunate ambiguity as to the physical significance of ‘‘ f”’ 
that, at first, made it seem to me impossible to accept and adopt Lacey’s 
equations. In this, I now think I was wrong, for the following reasons. 


(t) See p. 11. In Part I the basic dodge resorted to (in order to avoid the dilemma presented 
by our ignorance of the actual shape of the cross-section) was to assume an “ average ”’ cross-section, 
approximately intermediate, between the two extreme possibilities. See p. 19. 

(2) Lacey, 1930 ; Lacey, 1935. See Bibliography. 

(3) It is tolerably clear that they constitute the first major advance in our knowledge since 
Kennedy’s classic paper of 1895. As far as I personally am concerned, I owe the whole of what 
little I know about this most difficult subject, to a laborious and long continued effort to understand 
these two papers and the varied and’ voluminous criticism they evoked ; and to evaluate their 
importance, validity and shortcomings in the light of that criticism. I passed through the usual 
three stages of thinking (a) that they solved everything ; (5) that they solved nothing ; and (c) that 
they constitute a very real advance but that they still leave many problems unsolved. They leave 
standing, I fear, very few of the theoretical postulates that we owe to Gilbert’s two important papers 
(Gilbert, 1877 and 1914) and which geologists have almost universally accepted as valid. 

(4) i.e. physical significance. Cf. p. 17, note (4) infra. 

(5) Or at least one aspect of that load. 

(6) And is directly proportioned to /d ; when d=the mean diameter of the particles that make 
up the bed-load. 

(7) On this point, see especially the papers by Sir Claude Inglis in the Annual Reports of the 
Central Board of Irrigation, Poona. See also two papers by him which he tells me are shortly to 
appear in the Proc. Inst. Civ. Eng. 

(8) As regards the cross-section, Lacey maintains that a stream’s ultimate, graded cross-section 
would be a semi-ellipse. This is disputed by eminent authorities. See, for instance, Lane, 1937, 
p. 140. But both the cross-section assumed by Lacey and the suggested alternatives are—for our 
purposes—equivalent ; for they are all approximately equivalent to our own “ average cross-section ”’. 
See p.9. His assumptions as to cross-section, therefore, even if wrong, need not deter us from adopt- 
ing his equations. 
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It is certain that we do not yet know enough to calculate, even approxi- 
mately, the ultimate effect upon the hydraulic variables of a given change 
in the load or charge that a stream is called upon to transport”. We have, 
therefore, perforce, ignored the charge delivered for transport as an indepen- 
dent, external variable. In other words, we have tacitly assumed a constant 
charge, and this unfortunately we must continue to do until our knowledge 
is increased. 

We had better, therefore, make a merit of necessity and a iaowledee 
explicitly that our equations are based upon two assumptions : first, that the 
charge of detritus shall remain substantially constant ; second, that the mean 
erade of detritus shall remain substantially constant. 


This has the merit of intellectual honesty, and avoids self-deception ; 
and I have gradually come to see, first, that by this method we can get, straight 
from Lacey’s papers, and without becoming entangled in equations that 
contain the obnoxious “f”’, much (though not all) that we, as Geologists, 
require ; and second, that by no method, in the present state of our knowledge, 
can we get any more than this. 

Furthermore, a little reflection shows that of these two assumptions, the 
first at least is neither arbitrary nor improbable. ‘The primary determinants 
involved would seem to be (1) the amount of detritus available; (2) the 
gradient of the valley-side slopes ; (3) the denseness of the vegetable cover ; 
(4) to what extent there is “ storage ’’ of precipitation in the form of snow and 
ice ; (5) to what extent natural seepage is checked and run-off enhanced by 
frozen soil or subsoil; (6) the amount; and (7) the temporal distribution 
of the mean annual precipitation. 

Of these seven, the first five seem to be quite dehors *‘ Hydraulics ’’, even 
sensu latissimo ; while as regards (6) and (7) any change that tended to increase 
(or decrease) Q would also tend to increase (or decrease) the load delivered 
to the stream by its headwaters and valley-side slopes. ‘here would, therefore, 
be a tendency to maintain, at least roughly, a constant charge, despite 
substantial changes in rainfall, in Q,, and in the load®’ 

The first step towards the desired equations is this: An analysis of a 
mass of empirical data by the method of least squares led Lacey to conclude 


) Cf. p. 13 supra. By load I mean the total number of lbs. of detritus carried past a specified 
cross-section per second. ‘The charge=the load—Q. It is the amount (by weight) of detritus carried 
per cubic foot of water ; or, to be meticulously accurate, per cubic foot of (water plus detritus). Note 
that we are speaking of the charge that a stream is called upon to transport. Note, also, that the difficulty 
arises, primarily, with regard to the ultimate effects. If the stream is already graded and the charge 
is increased, the immediate effect on the variables is nil. Those who wish to retain the phrase “ an 
overloaded stream ”’ must tread warily. In these days “‘ self evident propositions ”’ are (rightly) at a 
discount. Our poor old friend the ‘‘ Excluded Middle ”’ has had a shocking time of it, but the Law 
of Contradiction still stands. And few propositions seem more certain than the assertion : a stream 
cannot transport more than it can transport. 


(2) Of the assumptions of constancy of charge and constancy of mean grade, as well as of the 
relations between mean precipitation, run-off, discharge, etc., I hope to say more on a subsequent 
occasion. 


16 


GRADED STREAMS 


that, in graded streams “of constant silt grade ’’, there exists between V and R 
the “ simple fundamental dynamic relationship”? : V,. « s/ Re 
@ 


oo Wg, ee TR a ae ue te ns (1) 


Similarly), empirical data also yielded, for graded streams ‘‘ of constant 


silt grade ’’, the relation : Sc oe 
G 


I 


: bi ef ee ge 
ite: Re ay a2 7 os me iB) 


Equation (1) expresses mathematically the well-known fact that, other 
things being equal, to increase the size of a stream) tends to increase its 
velocity. Likewise equation (2) expresses the easily verifiable fact that a 
small graded stream tends to have a steeper slope than a large one. Yet, 
mathematically, these two apparently innocuous equations immediately yield 
the apparent paradox that : 

I 


he et 
Vv = Re a be iy oo) = i 


If we interpret this equation as asserting (as it ‘‘ obviously ” appears to 
assert) that to increase the slope of a stream will tend to decrease its velocity, 
that proposition is (equally “ obviously’) nonsense. But if we merely 
transpose this equation, an operation devoid of any significance, mathematic- 
cally, we should write instead : | 


i 
So= Wye sins 36 Se °° 5 (4) 


And if we then interpret equation (4) physically, im conjunction with the 
primary equations from which we derived it, we see that it merely expresses 
a harmless tautology. It asserts, in effect, that since a relatively high V 
implies a relatively high R (equation (1)), therefore, if we are told that a 
graded stream has a relatively high V we may infer that it has a relatively 
high R, and may therefore infer also (equation. (2)) that its S will be relatively 
low™, 

For Lacey, the ‘fundamental variables’? are V, R and S Chand 


1) Lacey, 1935, p. 423. 

2) Lacey, 1935, p. 425. 

3) Whether graded or ungraded. 

4) I have stressed this point for two reasons. It can prove, as I know to my cost, a real difficulty 
to one seeking to understand.Lacey’s papers. It is also, or so it seems to me, a striking illustration of 
the important but insufficiently emphasized fact that in the interstices of ail mathematical symbolism 
there lurk indubitable gremlins. These are powerless to harm us so long as we stick to mathematics ; 
but can, and often do, cause us to talk the most shocking nonsense when we think we are, and appear 
to be, merely giving to the symbols of our final equation their ‘‘ obvious physical significance. 


(5) Lacey, 1935, p. 536. 


( 
( 
( 
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equations (1), (2) and (4) constitute the “‘ fundamental ” dynamic equations”, 
They express the interrelation that exists between the Sg, the Vg and the Rg of 
a graded stream carrying a constant charge of a constant mean grade of 
detritus. 

If we were told that a specified change in the regional rainfall would 
result in a quadrupled Rg, we should then know from equations (1) and (2) 
that the ultimate V< would be doubled and the ultimate S, halved: for 
equation (1) asserts that when the river has had time to re-grade itself, the 
dynamics of graded flow require that VF = R%* = 44 = 2; and equation (2) 


: | 
thatesc. = RH = a oe 


But the snag is that all three of these variables are dependent variables. 
S, is as much so. as Vg or Rg”. And unfortunately there is no one to volunteer 
such informatory titbits as to what the ultimate effect of a given change will 
be upon any one of these three interlocking dependent variables. Moreover, 
the ultimate shape of the channel has also changed sides and joined the 
category of dependent variables’. We are left, therefore, with but a single 
independent variable, Q'.. And we can make no use whatever of our three 
fundamental dynamic equations, until we have somehow derived an equation 
that will connect one of those interlocking variables with Q. 

Faced with this problem, Lacey turned for assistance to “the laws 
governing similarity’, a branch of mathematical physics that has been 
developed primarily in connection with the construction of scale models for 
experiments in miniature. 

His reasoning is as follows : We have already seen that in nature a small 
graded stream will have a steeper slope than a large one, if the two streams 


are in all other respects “‘ similar”. But if we think of the smaller stream as 
: (5) 

a natural model of the larger, and recall that S = [> we see that this entails 

a certain exaggeration of this ratio in the “ model”. In which case, to avoid 


“ planimetric distortion”) the laws of dynamic similarity require that all 
vertical dimensions in the ‘‘ model”? must be exaggerated, relative to its 
horizontal dimensions, to precisely the same extent. Thus, if the slope of 
the larger stream is 1 ft. per mile and that of the smaller 2 ft. per mile, all 
similar ratios in the latter, such as the steepness of the banks or the ratio of 
depth to width, must also be doubled. And Lacey tells us that Professor 


(1) Lacey, 1935, P- 437- 

(2)s5ee%p.07: ’ 

(3) Apart, of course, from the independent variable, “‘ charge delivered ”’, which we are forced to 
assume constant. 


(4) Lacey, 1935, Pp. 437- 
(5) Where h=height and 1=length. 


(6) Ibid., pp. 437 and 439. 
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Osborne Reynold, from his observations on natural streams, arrived, indepen- 
dently and empirically at the same conclusion™. 


eas R. 
But if this is in fact a natural ‘“‘ law ”), then the rest is easy : for psa 
precisely similar ratio toz and when we compare a stream regraded after a 


decrease in Q with that same stream before the decrease, we are in fact 


Rew he 
comparing two different streams. Hence : P ae See Dee 
fc) G 
I I 
But, from equation (4), S, « Weise Ng Vv. 
K pa 
- So Wn anion Vv. when K is a constant. 
PARE VE a RCV, 


But, from equation (1) V 0 4/R .. R o V? «. R = KV. 
Substituting in (5) : 


POE RG a hy cede (6) 


Substituting in (6) : 


Oc Ove oP is rod 
Ae SO ere Nye 


Cay ay c G me (52) 


OE TNE CNA TL uf 
But, by definition, of 10). aid eV Come ence: | 
s— ve" is 
i Wipe O eke Bh iM be he ne (8) 
Or, in Lacey’s own words?) : ‘‘ The author therefore concludes that in 


(1) Lacey, 1935, p. 437. Professor Reynold’s work, and the now famous “ Reynold’s Number ”’ 
are both too well known to require comment. But that he should have arrived inductively at the 
conclusion which Lacey reaches deductive ly is both interesting and important. 

(2) Lacey, 1935, p. 438. 

(3) Lacey, 1935, p. 438. 
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similar”) regime channels” of the same silt grades the wetted perimeters vary 
as the square roots of the discharges, and the velocities as the one-sixth power 
of the discharges. This relationship would appear to be a law.” 

To which it is only necessary to add the reminder that for the phrase 
‘‘ of the same silt grade ’’, it is probably necessary to substitute some wider 
phrase that will include both aspects of a river’s‘load—the charge as well as 
the mean grade of detritus. 

We have now all that we require: for equations (1) and (4) combined 
with equation (8) give us immediately equations expressing R& and Sé in 
terms of Q% ; while equations (6) and (8) give us Pé in terms of Qé. 
And since throughout this paper we are concerned solely with streams in 
which W>10 D, W = P (approx.), D = R (approx.), and D oc D®’. 

Our new, and I believe valuable, acquisitions, therefore boil down into 
the four amazingly simple equations : 


tek Ke 
VE= G 


( 
St = on (9) 
Bye ea tye <3 ae a, * | ( 
As to W we can even take the matter one step further. Lacey contends 
that in a graded stream, free to fashion a channel of a shape, depth and width 
of its own choosing, the ultimate P will be independent of “f” and will be 
determined solely by Q™. In this one case, therefofe, it is possible to give an 
equation that will yield an absolute value and not merely the value relative to 
the initial value which in our notation is symbolized by Wé&.  Lacey’s 


(1) If two streams (carrying equal charges of the same mean grade of detritus) have both, by 
natural processes, achieved a graded slope and generated an appropriate “ graded ”’ cross-section, 
they are “‘ similar’. Such streams do not exhaust the class: “similar streams”. Cf. Lacey, 1935, 
p. 438 et seg. But fortunately, if the proposition first stated is true, that is all that is requisite for our 
purposes and we need neither discuss generally nor attempt to define the concept “ similarity ”’. 

(2) Lacey, being an engineer, not a geologist, uses throughout both papers the expressions “ regime 
channel ’”’, “‘ regime flow’’, and the like. He often, however, uses instead the expression “ stable- 
rivers ” or “‘ stable channels ”’ (Lacey, 1930, p. 259; 1935, p- 421). He also speaks of “ the velocity 
at which equilibrium or stability is maintained ” (Lacey, 1930, p. 259) ; and (Lacey, 1936, p. 421) 
says: “‘ The problem of . . . regime flow . . . was one of silt transport by a channel that neither 
silted nor scoured, but swept forward an undiminished silt charge.” Again (Lacey, 1935, p. 426) 
he says: ‘*.There would appear to be a period . . . when such rivers, having ceased to scour, have 
not begun to deposit silt. ‘That period is synchronous with the peak of the flood, and the river is then 
momentarily in regime.” Such statements (of which many further examples occur throughout both 
papers) together with Kennedy’s original discussion (Kennedy, 1895) which Lacey expressly adopts 
as his starting point, make it quite certain, I think, that Lacey’s “‘ regime channel ”’ is synonymous 
with a geologist’s “‘ graded stream”. My sole reason for labouring the point is that my friend 
Sir Claude Inglis has at times suggested in discussion, that Lacey’s “‘ regime channel” is a more 
restricted concept than my “‘ graded stream ”’. But in what way more restricted, I have never been 
able to discover. 

isieoecenp.20. 

(4) Lacey, 1930, pp. 272-273; and 1935, p. 440. One of the premises on which Lacey’s equation 
rests is the hypothesis of a semi-elliptical cross-section. But the value yielded by the equation should 
still be of the right order of magnitude even if that hypothesis is untrue. Cf. p. 15, note 8 supra. 


é 
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equation may be written with sufficient accuracy for our purposes (‘in 
foot-second units’) : W = P = 8 QO}. 

Thus a graded stream that discharged, at bank-full stage, 10,000 cubic feet 
per second, should ultimately cut for itself a channel having a mean width 
of 3 X 100 = 300 feet. I leave it to my readers to examine for themselves 
how far the facts of nature tend to prove or disprove this simple and 
unambiguous deduction from Lacey’s fundamental hypotheses. For their 
information I might add that the Thames at Teddington, at bank-full stage, 
has a discharge of about this order of magnitude", 

In conclusion let us apply these new equations to the problem we con- 


sidered on page 10 of Part I. Since Q@=4 :— V#=4'=1-26; S* = i 
=0-79; Dé=4'=1-59; and Wé=4/4=2. So we now know—approxi- 
mately—what the graded reaches of the Thames must have been like towards 
the close of Atlantic times. We should have seen a river twice as wide and 
more than half again as deep as the present one, and flowing, on the average, 
about 25 per cent more swiftly ; yet a river the graded slope of which we 


should have found was only # as steep as that of the present Thames. 


(1) Thames Conservancy Statistics, pp. 20-21 ; 4,770 million gallons per day. 
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